Novel disordering mechanism in ferromagnetic systems with competing interactions 



(N ■ 

o : 

(N . 

>>: 



oo 



Juan Carlos Andresen, 1 Creighton K. Thomas, 2 'EI Helmut G. Katzgraber, 2 - 1 and Moshe Schechter 3 

1 Theoretische Physik, ETH Zurich, CH-8093 Zurich, Switzerland 
2 Department of Physics and Astronomy, Texas A&M University, College Station, Texas 77843-4242, USA 
3 Department of Physics, Ben Gurion University - Beer Sheva 84105, Israel 
(Dated: March 3, 2013) 

We study the interplay between ferromagnetic and spin-glass phases in magnetic systems such as 
LiHo x Yi_ x F4 theoretically, as well as using numerical simulations of a three-dimensional diluted long-range 
dipolar Ising model with competing interactions in the presence of a random field. Our results suggest the ex- 
istence of a novel disordering mechanism of the ferromagnetic phase due to the underlying spin-glass phase 
when a random field is applied. We numerically compute the zero-temperature phase boundary between the 
quasi-spin-glass and ferromagnetic phases as a function of the Ho concentration and random-field strength, and 
explain the peculiar linear dependence of the critical temperature on the strength of the random field found in 
recent experiments by Silevitch et al. [Nature 448, 567 (2007)] on the LiHo x Yi- x F4 compound. 
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Introduction. — The random-field Ising model (RFIM) 
plays a central role in the study of disordered systems and has 
been applied to problems across disciplines ranging from dis- 
ordered magnets to random pinning of polymers in two space 
dimensions, as well as water seepage in porous media. 

At and below the lower critical dimension di = 2, the fer- 
romagnetic (FM) phase is unstable to an infinitesimal random 
field US- At higher space dimensions the disordering of 
the ferromagnetic phase requires the random field strength to 
be of the order of the strength of the interaction between the 
spins. Yet, the effect of the random field on the phase tran- 
sition between the FM and paramagnetic (PM) phases — for 
systems with both short range and dipolar interactions — has 
been source of vast experimental and theoretical scrutiny over 
the last two decades; see Refs. and references therein. 

Over the last three decades the RFIM has been studied exper- 
imentally via dilute antiferromagnets in a field H, as both the 
random-field Ising model and the diluted antiferromagnetic 
Ising model in a field share the same universality class near 
criticality. More recently it has been shown that in anisotropic 
dipolar magnets the RFIM can be realized in the FM phase: 
By applying a transverse field to a dilute dipolar ferromag- 
net, such as LiHo x Yi_ x F4, one transforms the spatial disor- 
der to an effective random field in the longitudinal direction 
13-0]. This new possibility to study the RFIM in the FM 
phase opened the doors for advancing our understanding of 
the random field problem |Hoj|, as well as potential new ap- 
plications, such as the advent of tunable domain-wall pinning 
111 111 in magnetic materials. 

Indeed, Silevitch et al. recently studied the FM-to-PM 
phase transition in the presence of random fields in the 
LiHo x Yi_ x F 4 compound 112I1 . In particular, they found a 
peculiar behavior of the susceptibility as a function of the ap- 
plied transverse field: A singularity in the field dependence of 
the critical temperature at zero field, as well as a linear de- 
pendence of the transition temperature on the effective ran- 
dom field strength (which is linear in the applied field for 
small fields J3 0]). Similarly, in the molecular nanomagnet 



Mni2 — ac a strong suppression of the transition temperature 
as a function of the magnetic field was found II 1 3fl . However, 
in M1112 — ac it is unclear if for small fields the behavior is lin- 
ear as well and singular as a function of the transverse field. 

Both the LiHo x Yi_ x F4 and Mn.12 — ac systems, with small 
applied transverse magnetic fields, are well modeled by the 
RFIM. However, an important difference between the two sys- 
tems is that the random field in the former is a result of spatial 
dilution, dictating a broad distribution of the interaction terms 
in both magnitude and sign, whereas in the latter the random 
field is a consequence of small random tilts in the spin direc- 
tion. Therefore, it is unclear which "physical ingredients" are 
required such that the transition temperature depends linearly 
on the random field strength, and, in particular, which under- 
lying mechanism dictates this behavior. 

Here we study the interplay between FM and spin-glass 
(SG) phases in a dipolar Ising model with competing inter- 
actions in the presence of a random field. We find a novel 
disordering mechanism of the FM phase when a random field 
is applied and the system is in close proximity (e.g., via dilu- 
tion) to an underlying SG phase. This disordering mechanism 
lies somewhat between the Imry-Ma and standard disorder- 
ing mechanisms: The disordering of the FM phase occurs at 
a finite field, which is considerably smaller than the typical 
spin-spin interaction strength. At zero temperature we predict 
the existence of a FM-to-quasi-SG phase transition and nu- 
merically determine for LiHo x Yi_ x F4 the phase boundary as 
a function of the Ho concentration x and random field strength 
h. At finite temperature our theory is compatible with exper- 
iments at various Ho concentrations lfl2ll . suggesting that it is 
indeed the existence of competing interactions and the prox- 
imity to the SG phase that dictate the peculiar dependence of 
the critical temperature T c on the random field strength h. 

Theoretical analysis. — We first study the dilute 
LiHo x Yi_ x F4 system at zero temperature. For dilu- 
tions x > x c the system is FM, whereas for xq < x < x c the 
system is a spin glass (SG). Below we show numerically that 
x c rj 0.3. To date, it is unclear if xq > HQ. For x ~ x c 
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we define the energy per spin of the lowest ferromagnetic 
state of the system as /fm(^); the lowest energy of the spin- 
glass state as fsoix). Note that /pm(^) is the ground-state 
energy of the FM phase when x > x c and fsc(x) represents 
the ground-state energy of the SG phase for xo < x < x c . At 
x = x c /fm(i c ) = /sg(^c)> an d for x rss x c , to first order 
in x - x c , fsa(x) - Ifm(x) = a(x - x c ) + .... Consider 
the FM phase for x > x c and apply a random field of average 
strength h. For small random-field strength, the FM state in 
three dimensions cannot gain energy from the field, because 
domain flips are not energetically favorable. However, 
for spin glasses the lower critical (Imry-Ma) dimension is 
infinity IU6I1 (no order in a field). In particular, in three 
space dimensions the energy of the system can be lowered by 
flipping domains, creating a quasi-spin-glass (QSG) phase 
with a finite correlation length. Thus, for x close enough to x c 
the energy of the SG state will be lower than the energy of the 
FM state at a finite random field, which is still considerably 
smaller than the typical interaction strength between the 
spins. More generally, any three-dimensional Ising system 
with competing interactions having at zero random field a FM 
ground state and a SG state in a somewhat higher energy, will 
have a phase transition to the QSG phase as a function of the 
random field strength. The disordering of the ferromagnetic 
phase therefore occurs at a finite field whose magnitude de- 
pends on the proximity to the SG phase and, in particular, can 
be much smaller than the characteristic spin-spin interaction 
strength. Because in systems like LiHo x Yi_ x F4 the effective 
random fields are a result of quantum fluctuations fl7ll . 
this phase transition is a particular case of a quantum phase 
transition where the quantum fluctuations of the individual 
spins are small involving only the spin's ground and first 
excited states 11811 . but where the collective effect of all spins 
is strong enough to drive the transition. 

The value of the critical random field can be estimated us- 
ing the following Hamiltonian that describes a short-range 
Edwards-Anderson Ising spin glass lfl9ll in a random field: 
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Hea — — ^ JijSiSj + ^ hiSi. 

(ij) » 



(1) 



In Eq. ([]]) Jij are nearest-neighbor Gaussian random bonds 
between the spins Si with zero mean and standard deviation J, 
and hi are Gaussian random fields of average strength h J20I1 . 
The SG ground state is unstable to an infinitesimal random 
field, creating domains of typical size (J ' /h) 1 ^ z / 2 ~ s \ with 
9 « 0.2 the stiffness exponent in three dimensions 112 111 . The 



energy reduction per spin resulting from the random field is 
thus of order /(/i) = h( J /h)^' 3 /^/^/ 2 -^. One can show 
that the total energy reduction per spin is of the same order, 
because the energy cost to flip domains is much smaller. The 
critical field h c (x) can be computed from f(h = h c ) = fsG~ 
/fm = ol(x — x c ) to obtain 



h c (x) = a'(J)(x - x c )^^ , 
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(2) 
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where a' (J) = a 3 - e Jz=v, see Fig. [1] For h > h c (x) 
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FIG. 1: (Color online) Comparison of the zero-temperature numer- 
ical and analytical [Eq. ifJJ] phase diagrams for the diluted dipolar 
Ising model [Eq. ©J in the h-x plane. Using 8 ~ 0.2 I2H1 we obtain 
h c (x) ~ (x — x c ) 464 . The analytical prediction, Eq. ((2), agrees 
very well with the numerical data with a' the only fitting parameter. 



the system consists of finite-size domains within which short- 
range glassy order persists. The size of these domains de- 
creases with increasing field, where at h « J the system re- 
sembles a simple paramagnet. As x — > x c the disordering 
field approaches zero. For large x — x c and h « J there is a 
crossover to the standard behavior where the disordering is a 
result of single-spin energy minimization, i.e., the intermedi- 
ate QSG regime disappears. 

When both T and x x c tire small but finite and the ran- 
dom field strength increases, the system changes first from the 
FM phase to the QSG phase, where the system is still dynam- 
ically frozen. Increasing the random field further, the QSG 
domain size decreases. The crossover to a disordered para- 
magnet is indicated by a peak in the magnetic susceptibility, 
which is broad and of finite magnitude, its width dictated by 
the typical domain size upon melting f7[]. This scenario has 
been observed experimentally; see Fig. 2 of Ref. Ill 21 . As the 
temperature is increased, the crossover to the PM regime oc- 
curs at smaller random fields, dictating larger typical domains 
and a higher peak for the magnetic susceptibility Ip7 , 22 1 . 

Let us denote by Tsg the ordering temperature at zero ran- 
dom field of the underlying SG phase. For T > TgG the 
domain size at the crossover to the QSG phase is limited by 
£fm oc (T — Tsg) -1 '. where v is the critical exponent of the 
correlation length. When the crossover field is small enough 
for £fm to be reached, further increase in temperature leads 
to a smaller critical field h c . However, the magnitude of the 
susceptibility peak saturates; see also Fig. 2 of Ref. |[l2"ll . 

To quantify the above observations and analyze the depen- 
dence of T c on the effective random field, let us denote the 
lowest free energies of the FM phase below the critical FM 
temperature Tfm and a competing QSG phase, in which SG 
correlations exist within domains of size £, as Fpm{x,T) 
and Fqsq(x,T), respectively. Because the entropy of the 
QSG phase is dominated by the regions at the boundaries be- 
tween domains lfl6ll . the main effect of the random fields is 
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to lower the energy of the QSG phase. Thus, Fqsg(%, T) 
F FM (x, T) = —A(T - T FM ) + B(x - x c ) - h/^ G . He 
we used the fact that for h = Fqsq(x c , T c ) = Ffm(x c , T ( 
We thus obtain 



T c (h) - T c (0) 



1 

A 



B(x-x c ) -h/i 



3/2 



(3) 



Equation ([3]i assumes that the FM phase is disordered by a 
QSG phase, in which domains of size £ exist. Thus, it gives 
an upper bound for the value of T c {h). In fact, as we ex- 
pect T c (h) to be a monotonous function, we expect the above 
expression for T c (h) to be valid at ft, > h* , where h* — 
B{x-x c )f/ 2 is defined by the condition T c (h) = T c (0). Our 
theoretical interpretation of the experimental data in Ref . lll2"ll 
is in good agreement with our numerical results presented be- 
low for the LiHo x Yi_ x F4 system at different dilutions and 
magnetic fields. Both our theoretical and numerical results are 
consistent with T c (h) being linear in a regime where h <^ J, 
with deviation from linearity as h — > 0. Note that in Ref. 11211 
T c (h) is linear down to the lowest random fields studied if one 
defines the critical temperature by the asymptotic behavior of 
the susceptibility at high temperatures. However, if T c (h) is 
defined by the peak position of the susceptibility, deviations 
from linearity are observed at low fields 023J. 

Numerical details. — Previous work showed that 
LiHo x Yi_ x F4 at low temperatures and in the presence of an 
external transverse magnetic field is described by an effective 
Hamiltonian consisting of three terms Jg, 24]: a long range 
term [dipolar, first term in Eq. (0}], an antiferromagnetic 
short range term [second term in Eq. (0}], and a longitudinal 
random field term [third term in Eq. (@)]. We thus study 
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Here ej = {0, 1} is the occupation of the magnetic Ho 3+ ions 
on a tetragonal lattice (lattice constants a = b = 5.175A and 
c = 10. 75 A) with four ions per unit cell SHil. Si G {±1} 
are Ising spins, hi represents the random fields drawn from a 
Gaussian distribution with zero mean and standard deviation 
h, where h is the strength of the applied field in [K]. The 
magnetostatic dipolar coupling Jy between two Ho 3+ ions is 
given by: = D(r?- - 3z?-)/rf •, where = |r, - rj\, r, 



is the position of the i-th Ho 3+ ion and 
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the component parallel to the easy-axis. The dipolar constant 
is D/a 3 = 0.214K (27ll and the antiferro mag netic nearest- 
neighbor exchange is set to J ox = 0.12K 125l 12811 . We use 
periodic boundary conditions. The long-range interactions are 
taken into account via Ewald sums 12a l29ll . At zero field and 
no dilution we find T c = 1.5316(2)K, in excellent agreement 
with experimental results where T c = 1.530(5)K 1 30 ] . 

For the zero-temperature simulations [Fig. Q] we use 
(heuristic) jaded extremal optimization Oil 13211 . For each 
run, we use exponents t = 1.6, 1.8 and 2 with an aging pa- 
rameter r = 0.05 for at least 2 26 steps. Empirically, ground 
states are found with high confidence for L < 10 and h = 0, 



and L < 8 with small h ^ 0. Because the ground states 
of this model tend to be unique, we use the Binder ratio to 
identify the phase boundary. The Binder ratio is defined by 

g = (1/2) (3 - [m 4 ] av / [m 2 ] 2 J, where m = (1/^)^5, 
(N = 4xL 3 the number of spins) is the magnetization of 
the system and [■ • •] represents an average over the disor- 
der. The Binder ratio is a dimensionless function and scales 
as g ~ G\L X I V (x — x c )], allowing for the extraction of x c 
and v for a fixed value of h. Simulation parameters are listed 
in the supplementary material section, Table U 

At finite temperatures we use the parallel tempering Monte 
Carlo method 113 3H . Simulation parameters are listed in Tables 
[II] [TTT| and [IV] in the supplementary material section. To deter- 
mine the finite-temperature transitions for a given value of h 
and x we measure the two-point correlation length II34|| 



1 / l(m 2 (0)) T ], 

2sin(fc min /2)\/ [(m 2 (k min )) g 



- 1 



where 



1 N 

i=i 



kR, 



(5) 



(6) 



Here (■ • -) T represents a thermal average, and R; is the spa- 
tial location of the spin Si. k m i n = (2vr/i,0, 0) represents 
the smallest nonzero wave vector. Near the transition is 
dimensionless and expected to scale as £l / L ~ X[L 1 / iy (T — 
Tc)]. Because at the transition temperature, T = T c , the argu- 
ment of the scaling function is zero (up to scaling corrections) 
and hence independent of L, we expect lines of different sys- 
tem sizes to cross at this point [Fig. Ufa)]. If however the lines 
do not meet, we know that no transition occurs in the stud- 
ied temperature range [Fig. 0c)]. To determine the critical 
temperature T c (h) we perform a finite-size scaling analysis of 
the data [see Fig. Efb)] . Using a Levenberg-Marquardt mini- 
mization combined with a bootstrap analysis as described in 
Ref. 113511 allows us to determine the optimal values of the crit- 
ical parameters with a statistical error bar, see Table IV1 

Figure [T] shows the h-x phase diagram of the 
LiHo x Yi_ x F4 system at zero temperature. We find 
very good agreement with Eq. (|2j, using 8 w 0.2 12 ill , i.e 
h c (x) ~ (x — x c )° , and a' as a fitting parameter. Figure|3] 
shows finite-temperature data for different Ho-concentrations 
x. Panel Oa) shows the critical temperature T c as a function 
of the random field h for x — 0.32, i.e., x — x c = 0.02 
small. Our results at finite T corroborate our theoretical 
model where the FM phase disorders at h w 0.045(5), a 
value slightly larger than found with the T = simulations. 
We also find that T c (h) at low fields is well described by a 
linear behavior with possible deviations from linearity at very 
small fields, as suggested by our analytical considerations 
above. Both the disordering of the FM phase at small fields 
and the linearity of T c {h) seem to persists up to a; = 0.44 
[Fig-Ob)], the dilution used in the experiments of Silevitch et 
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FIG. 2: (Color online) Finite-temperature data for the finite-size correlation length £l/L as a function of temperature T and x = 0.32. (a) 
Data for h = 0. There is a clear crossing for T c — 0.340(3) (2^1 for different system sizes L. (b) Scaling analysis of the data for h = 0. The 
solid line corresponds to the optimal scaling function (polynomial approximation), (c) Data for h — 0.06. There is no sign of a transition. 
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FIG. 3: (Color online) Critical temperature T c 12811 as a function of the random field strength h for different dilutions x. (a) x = 0.32. The 
critical temperature is linear in h. Moreover, above h « 0.045(5) there is no sign of a phase transition, (b) Similar behavior as for x = 0.32 
is found for x = 0.44 (dilution used in the experiments of Ref. fl2ll ). (c) x = 0.65. At large dilutions T c (h) is quadratic in h. 



al. II 1 211 - For x = 0.65, i.e., far from the SG phase, [Fig. [3jc)] 
the behavior of T c (h) changes drastically to a clear quadratic 
dependence, suggesting a standard FM-PM transition at this 
concentration. Estimates of the critical parameters are listed 
in the supplementary material, Table IV1 

Note that our zero and finite-temperature analyses for zero 
random field suggest a temperature-independent critical con- 
centration separating the FM and SG phases, i.e., the phase 
diagram for the LiHo x Yi_ x F4 system as a function of con- 
centration and temperature does not contain a reentrant SG 
phase at zero field, in contrast to previous suggestions in the 
literature Q6I1 . At the same time, our results at finite random 
fields at zero and finite temperatures for x = 0.44 and 0.65 
suggest that there is a range of random fields where the sys- 
tem shows reentrance to a SG phase at low temperatures, a 
common effect found in glassy spin systems 13711 . 



Conclusions. — We propose a novel disordering mecha- 
nism for three-dimensional ferromagnets with competing in- 
teractions and an underlying spin-glass phase. Furthermore, 
we explain various aspects of the experimental results of Sile- 
vitch et al. jl2j| , including the peculiar linear dependence of 
the critical temperature on the applied transverse field. We 
further find that at smaller concentrations (x = 0.32, close 
to the spin-glass phase) the reduction of T c with the random 
field becomes more pronounced. Our results strongly support 
the notion that it is the interplay between the competing inter- 



actions and the induced effective random field that dictate the 
behavior of the LiHo x Yi„ x F4 ferromagnet at low concentra- 
tions. It would be interesting to study the generality of the 
above arguments with other anisotropic dipolar systems. 
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TABLE II: Simulation parameters at finite temperature and x = 
0.32 for different fields h and system sizes L. The equilibra- 
tion/measurement times are 2 b Monte Carlo sweeps. T mln [T max ] 
is the lowest [highest] temperature used and Nt is the number of 
temperatures. iV sa is the number of disorder realizations. 
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500 


15 


750 


n ^9 


020 


p. 


lO 


u.uou 


u.ouu 


OU 


IOUU 


n 


u.uzu 


o 
o 




n n^n 
u.uou 


u.ouu 


OU 


1 nnn 

1UUU 


0.32 


0.020 


10 


16 


0.245 


0.500 


15 


1000 


32 


020 


1 2 


10 


245 


500 


15 


000 


32 


020 


14 


17 


274 


500 


13 


7500 


0.32 


0.030 


Q 


14 


0.226 


0.450 


14 


3000 


32 


030 




17 


212 


450 


15 


2000 


32 


030 


10 


17 


226 


450 


14 


2000 


32 


030 


1 2 


19 


226 


450 


14 


000 


0.32 


0.035 





13 


0.218 


0.450 


21 ) 


3000 


0.32 


0.035 


8 


14 


0.218 


0.450 


20 


2000 


0.32 


0.035 


10 


15 


0.218 


0.450 


2(1 


2000 


0.32 


0.035 


12 


18 


0.218 


0.450 


20 


650 


0.32 


0.040 


6 


14 


0.218 


0.450 


20 


3000 


0.32 


0.040 


8 


16 


0.218 


0.450 


211 


2800 


0.32 


0.040 


10 


17 


0.218 


0.450 


20 


2000 


0.32 


0.040 


12 


17 


0.218 


0.450 


14 


1000 


0.32 


0.050 


6 


13 


0.100 


0.650 


35 


4500 


0.32 


0.050 


8 


16 


0.100 


0.650 


35 


3000 


0.32 


0.050 


10 


18 


0.100 


0.650 


35 


1000 


0.32 


0.060 


6 


15 


0.100 


0.650 


35 


4000 


0.32 


0.060 


8 


16 


0.150 


0.500 


20 


1000 


0.32 


0.060 


10 


IS 


0.150 


0.500 


20 


1000 


0.32 


0.060 


12 


17 


0.150 


0.500 


2(1 


1000 



Supplementary Material: Andresen et al. 

TABLE I: Simulation parameters for T — 0: System of size L = 6, 
8, and 10, field h and dilution x are studied. x m i n [a; max ] is the 
smallest [largest] concentration studied and Ax is the step size be- 
tween measurements. iV sa is the number of disorder realizations. 



h iCmin ^max As JV aa 

OOOO 0280 O350 O010 5000 

0.025 0.275 0.400 0.025 3000 

0.050 0.300 0.400 0.025 3000 

0.075 0.300 0.600 0.050 1500 

0.100 0.400 0.800 0.100 1500 



6 



TABLE III: Simulation parameters at finite temperature and x = 



0.44. For details see Table|Il] 


X 


h 


L 


6 


T min 


T mi „ 


N T 


iV sa 


0.44 


0.000 


6 


10 


0.500 


1.000 


30 


1000 


0.44 


0.000 


8 


12 


0.500 


1.000 


30 


1300 


0.44 


0.000 


10 


13 


0.500 


1.000 


30 


2500 


0.44 


0.000 


12 


14 


0.500 


1.000 


30 


550 


0.44 


0.000 


14 


14 


0.560 


1.000 


47 


650 


0.44 


0.020 


6 


10 


0.525 


0.800 


30 


2000 


0.44 


0.020 


8 


12 


0.525 


0.800 


30 


2000 


0.44 


0.020 


10 


13 


0.525 


0.800 


30 


1000 


0.44 


0.020 


12 


14 


0.525 


0.800 


30 


1000 


0.44 


0.020 


14 


14 


0.550 


0.800 


30 


900 


0.44 


0.040 


6 


11 


0.525 


0.750 


22 


2000 


0.44 


0.040 


8 


12 


0.525 


0.750 


22 


2000 


0.44 


0.040 


10 


14 


0.540 


0.750 


20 


1000 


0.44 


0.040 


12 


15 


0.550 


0.750 


20 


1000 


0.44 


0.060 


6 


11 


0.500 


0.750 


25 


2000 


0.44 


0.060 


8 


13 


0.500 


0.750 


25 


2000 


0.44 


0.060 


10 


13 


0.519 


0.730 


19 


2000 


0.44 


0.060 


12 


14 


0.519 


0.730 


19 


1300 


0.44 


0.080 


6 


12 


0.475 


0.725 


20 


2000 


0.44 


0.080 


8 


13 


0.475 


0.725 


20 


2000 


0.44 


0.080 


10 


13 


0.500 


0.725 


20 


1200 


0.44 


0.080 


12 


14 


0.509 


0.725 


19 


1300 


0.44 


0.100 


6 


13 


0.445 


0.725 


23 


2000 


0.44 


0.100 


8 


14 


0.445 


0.725 


23 


2300 


0.44 


0.100 


10 


13 


0.445 


0.725 


23 


2200 


0.44 


0.100 


12 


15 


0.475 


0.725 


20 


880 


0.44 


0.120 


6 


11 


0.445 


0.725 


23 


2500 


0.44 


0.120 


8 


12 


0.445 


0.725 


23 


2200 


0.44 


0.120 


10 


13 


0.445 


0.725 


25 


1000 


0.44 


0.120 


12 


14 


0.445 


0.725 


25 


1000 



TABLE V: Critical parameters estimated using a finite-size scaling 
technique: For each concentration x and field strength h we compute 
the critical temperature T c and critical exponent v. 



X 


h 


7c 


V 


0.32 


0.000 


0.340(3) 


0.83(3) 


0.32 


0.005 


0.335(3) 


0.89(3) 


0.32 


0.010 


0.340(3) 


0.81(3) 


0.32 


0.020 


0.314(5) 


1.04(7) 


0.32 


0.030 


0.281(13) 


1.23(17) 


0.32 


0.035 


0.269(14) 


1.31(19) 


0.32 


0.040 


0.247(14) 


1.41(15) 


0.44 


0.000 


0.584(1) 


0.75(1) 


0.44 


0.020 


0.581(1) 


0.70(1) 


0.44 


0.040 


0.563(3) 


0.77(2) 


0.44 


0.060 


0.548(5) 


0.86(3) 


0.44 


0.080 


0.522(5) 


0.99(4) 


0.44 


0.100 


0.506(5) 


1.01(4) 


0.44 


0.120 


0.466(9) 


1.39(12) 


0.65 


0.000 


0.9597(8) 


0.79(2) 


0.65 


0.050 


0.9531(10) 


0.76(2) 


0.65 


0.010 


0.9264(14) 


0.84(3) 


0.65 


0.015 


0.8832(21) 


0.91(4) 


0.65 


0.020 


0.8312(41) 


1.06(10) 


0.65 


0.030 


0.6905(113) 


1.11(12) 



TABLE IV: Simulation parameters at finite temperature and x = 
0.65. For details see Table HT1 



X 


h 


L 


6 


T min 


T min 


N T 


iV sa 


0.65 


0.000 


6 


10 


0.500 


1.400 


20 


1000 


0.65 


0.000 


8 


12 


0.500 


1.400 


20 


500 


0.65 


0.000 


10 


12 


0.500 


1.400 


20 


450 


0.65 


0.000 


12 


11 


0.500 


1.400 


20 


500 


0.65 


0.000 


14 


10 


0.850 


1.400 


15 


490 


0.65 


0.000 


16 


11 


0.850 


1.400 


15 


470 


0.65 


0.050 


6 


8 


0.697 


1.400 


15 


1000 


0.65 


0.050 


8 


10 


0.697 


1.400 


15 


500 


0.65 


0.050 


10 


9 


0.800 


1.400 


20 


500 


0.65 


0.050 


12 


11 


0.697 


1.400 


15 


300 


0.65 


0.050 


14 


10 


0.900 


1.400 


15 


500 


0.65 


0.050 


16 


11 


0.920 


1.400 


16 


280 


0.65 


0.010 


6 


8 


0.820 


1.400 


20 


750 


0.65 


0.010 


8 


9 


0.820 


1.400 


20 


500 


0.65 


0.010 


10 


10 


0.820 


1.400 


20 


500 


0.65 


0.010 


12 


11 


0.820 


1.400 


20 


800 


0.65 


0.010 


14 


12 


0.820 


1.400 


20 


380 


0.65 


0.015 


6 


8 


0.820 


1.400 


20 


1000 


0.65 


0.015 


8 


10 


0.820 


1.400 


20 


500 


0.65 


0.015 


10 


11 


0.820 


1.400 


20 


500 


0.65 


0.015 


12 


12 


0.820 


1.400 


20 


500 


0.65 


0.015 


14 


13 


0.820 


1.400 


20 


500 


0.65 


0.020 


6 


10 


0.661 


1.400 


15 


500 


0.65 


0.020 


8 


11 


0.661 


1.400 


15 


400 


0.65 


0.020 


10 


14 


0.756 


1.400 


13 


1300 


0.65 


0.020 


12 


15 


0.756 


1.400 


13 


1000 


0.65 


0.020 


14 


16 


0.756 


1.400 


13 


990 


0.65 


0.030 


6 


10 


0.600 


1.400 


20 


750 


0.65 


0.030 


8 


15 


0.600 


1.400 


20 


1000 


0.65 


0.030 


10 


17 


0.600 


1.400 


20 


550 


0.65 


0.030 


12 


19 


0.600 


1.400 


20 


530 



